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ABSTRACT 

This  paper  presents  an  analytical  simulation  of  the  dynamics  of  a  gear 
transmission  system  compared  to  experimental  results  from  a  gear  noise  test 
rig  at  the  NASA  Lewis  Research  Center.  The  analytical  procedure  developed 
couples  the  dynamic  behaviors  of  the  rotor-bearing-gear  system  with  the 
response  of  the  gearbox  structure.  The  modal  synthesis  method  is  used  in 
solving  the  overall  dynamics  of  the  system.  Locally  each  rotor-gear  stage  is 
modeled  as  an  individual  rotor-bearing  system  using  the  matrix  transfer  tech¬ 
nique.  The  dynamics  of  each  individual  rotor  are  coupled  with  other  rotor 
stages  through  the  nonlinear  gear  mesh  forces  and  with  the  gearbox  structure 
through  bearing  support  systems.  The  modal  characteristics  of  the  gearbox 
structure  are  evaluated  using  the  finite  element  procedure.  A  variable  time 
stepping  integration  scheme  is  used  to  calculate  the  overall  time  transient 
behavior  of  the  system  in  modal  coordinates.  The  global  dynamic  behavior  of 
the  system  is  expressed  in  a  generalized  coordinate  system.  Transient  and 
steady  state  vibrations  of  the  gearbox  system  are  presented  in  the  time  and 
frequency  domains.  The  vibration  characteristics  of  a  simple  single  mesh  gear 
noise  test  rig  is  modeled.  The  numerical  simulations  are  compared  to  experi¬ 
mental  data  measured  under  typical  operating  conditions.  The  comparison  of 
system  natural  frequencies,  peak  vibration  amplitudes,  and  gear  mesh  frequen¬ 
cies  are  generally  in  good  agreement. 
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NOMENCLATURE 

modal  function  of  the  ith  mode  in  x-direction 

modal  function  of  the  ith  mode  in  8-direction 

modal  function  of  the  ith  mode  in  y-direction 

gearbox  damping  matrices 
torsional  damping  matrix 

bearing  direct  and  cross-coupling  damping  matrices 
bearing  excitation  forces 
gear  mesh  torque 

gear  mesh  force  in  x-  and  y-directions 
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external  excitation  moment 
external  excitation  forces 
gryoscopic-angular  acceleration  matrix 
gyroscopic-angular  rotation  matrix 
identity  matrix 

rotational  mass  moment  of  inertia  matrix 
average  stiffness  matrix 
gearbox  stiffness  matrix 

compensation  matrices  in  x-  and  y-direction 
shaft  stiffness  matrix 
torsional  stiffness  matrix 

gear  mesh  stiffness  between  ith  and  kth  rotor 
bearing  direct  and  cross-coupling  stiffness  matrix 
mass-inertia  matrix  of  rotor 
mass-inertia  matrix  of  gearbox 

the  number  of  modes  used  to  define  each  motion 
radius  of  gear  in  the  ith  rotor 
gear  generated  torque 

generalized  motion  in  x-  and  y-directions 
gearbox  motion  in  x-,  y-  and  z-directions 
gearbox  motion  at  bearing  supports 

gear  displacements  in  x-  and  y-directions  of  the  ith 
rotor 

gear  forces  in  x-  and  y-directions 

motion  of  rotor  at  bearing  support 

angle  of  tooth  mesh  between  kth  and  ith  rotor 

rotational  displacement  of  the  ith  rotor  at  the  gear 

rotational  displacement  vector  of  the  ith  rotor 
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lateral  and  torsional  eigenvalue  diagonal  matrices 

lateral  and  torsional  orthonormal  eigenvector  matrices 
of  the  kth  rotor 

jth  orthonormal  mode  of  kth  rotor  at  1th  station 


INTRODUCTION 

The  study  of  gearbox  vibration  is  an  important  area  of  research  among  the 
engineering  community.  With  increased  speed  and  torque  requirements  and 
reduced  allowable  vibration  levels,  a  critical  need  exists  to  develop  an 
accurate  analytical  simulation  of  the  dynamics  of  gear  transmission  systems. 
Recently,  a  number  of  analytical  modeling  and  simulation  procedures  have  been 
developed  to  predict  the  dynamics  of  multistage  gear  transmission  systems. 

Some  experimental  work  has  also  been  performed  to  improve  the  understanding  of 
gear  stress  and  vibration.  There  exits  a  need  to  correlate  analytical  and 
experimental  efforts  so  that  the  analytical  methods  can  be  refined  and 
verified. 

There  is  a  wealth  of  literature  on  vibration  analysis  and  simulation  of 
gear  transmission  systems.  August  (1986)  simulated  the  vibration  characteris¬ 
tics  of  a  three-planet  transmission  system.  Pike  (1987)  and  later  Choy 
(1988b)  used  the  gear  relationships  introduced  by  Cornell  (1981)  to  study  the 
dynamics  of  coupled  gear  systems.  Mitchell  (1987)  applied  and  matrix  transfer 
method  to  simulate  gear  vibrations.  Ozguven  and  Houser  (1988)  and  Kahraman 
(1990)  used  a  finite  element  model  to  predict  the  dynamics  of  a  multigear  mesh 
system.  Choy  (1991)  applied  a  modal  synthesis  technique  in  conjunction  with 
finite  element  and  the  matrix  transfer  methods  in  both  time  and  frequency 
domains  to  calculate  the  transient  rotor  and  casing  motions. 

There  are  few  studies  correlating  analytical  predictions  and  experimental 
results.  Some  experimental  correlation  work  by  Lim  (1990)  compared  analytical 
predictions  of  vibration  of  a  gear  housing  and  mounting  with  measurements. 

This  paper  reviews  the  development  of  a  global  dynamic  model  of  a  gear 
train  system,  and  presents  results  of  a  comparison  made  between  the  analytical 
predictions  of  the  model  and  experimental  results  from  a  gear  test  rig.  A 
combined  approach  using  the  modal  synthesis  and  finite  element  method  for  ana¬ 
lyzing  the  dynamics  of  gear  systems  was  developed.  The  method  includes 
dynamic  coupling  between  the  housing  and  the  gear-rotor  system.  The  single 
mesh  gear  noise  test  rig  at  NASA  Lewis  was  used  as  an  example  in  this  analy¬ 
sis.  The  gearbox  modal  characteristics  were  evaluated  by  solving  a  finite 
element  model.  These  modal  characteristics  were  compared  to  those  measured  in 
the  experimental  study  to  insure  accurate  modeling  of  the  experimental  rotor- 
gear-housing  system.  Transient  vibrations  of  the  rotor-gear  stages  and  the 
gearbox  housing  were  analytically  evaluated  by  the  modal  method.  A  numerical 
FFT  (Fast  Fourier  Transform)  algorithm  was  used  to  examine  vibration  in  the 
frequency  domain.  Frequency  spectra  of  the  predicted  gearbox  vibration  were 
compared  to  those  measured  in  the  experimental  study.  The  following  discus¬ 
sion  and  conclusions  are  drawn  from  these  comparisons  between  the  experimental 
results  and  analytical  predictions. 


□  □ 


ANALYTICAL  PROCEDURE 


The  equations  of  motion  for  a  single  mesh  multimass  rotor-bearing-gear 
shaft  system,  with  the  coupling  effects  of  gearbox  vibrations  and  the  rotor 
inertia-gyroscopic  effects,  can  be  written  in  matrix  form  for  the  ith  rotor 
(Choy  1987,  1989)  for  the  X-Z  plane  as: 
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Here  F  and  F  are  force  excitations  from  the  effects  of  mass  imbalance 
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and  shaft  residual  bow  in  the  X-  and  Y-directions .  F  and  F  are  the  X 

Gx  Gy 

and  Y  gear  mesh  forces  induced  from  gear  teeth  interactions.  Tne  bearing 
forces  are  evaluated  from  the  relative  motion  between  the  rotor  {X},  {Y}  and 
the  gearbox  {Xb>,  {Yfe}  at  the  bearing  locations  (Choy  1987).  The  mass-inertia 
and  gyroscopic  effects  are  incorporated  in  the  mass  matrix  [M]  and  the  gyro¬ 
scopic  matrices  [Gv]  and  [Gft].  The  coupled  torsional  equations  of  motion  for 
the  single  gear-rotor  system  can  be  written  as: 

[4(4  >[<4(4  ♦[KT]i(0}i  =  (FT(t)}i  +  (F^(t)l  (3) 

In  Eq.  (3),  {FT(t)>  represents  the  externally  applied  torque  and  {Fct(t)> 
represents  the  gear  mesh  induced  moment.  Note  that  Eqs.  (1)  to  (3)  repeat  for 
each  single  gear-rotor  system.  The  gear  mesh  forces  couple  the  force  equa¬ 
tions  of  each  gear-rotor  system  to  each  other  as  well  as  the  torsional  equa¬ 
tions  to  the  lateral  equations  (Choy,  1989;  Cornell,  1981;  David,  1987  and 
1988).  In  addition,  there  are  equations  of  motion  for  the  gearbox  which 
couple  the  various  gear-shaft-systems  through  the  bearing  supports.  The  gear¬ 
box  equations  can  be  written  as: 

X-equation 

MW  *  MM  *  MK*-  -  M  *  Ml*-  -  M 4  MINI 

(4) 
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Y-equation 
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N(M  -  M(M  ♦  [=,»]{*»  -  M * [c,y]{4»  - *.}  *  Ml1.) 


(5) 


and  Z -equation 


[Mb](2b}  +  [Cbz]{2b}  +  [Kbz]{Zb}  -  *b,<*> 


where  F  (t)  is  an  excitation  function  due  to  external  forces  in  the  axial 
bz 

direction.  Since  the  bearing  is  assumed  to  be  uncoupled  in  the  Z-direction, 
Eq.  (6)  can  be  solved  independently  without  considering  shaft  motion. 

The  torsional  and  lateral  vibrations  of  a  single  individual  rotor  and  the 
dynamic  relationships  between  each  gear-rotor-system  are  coupled  through  the 
nonlinear  interactions  in  the  gear  mesh.  The  gear  mesh  stiffness  varies  in  a 
periodic  nonlinear  pattern  with  each  tooth  pass  engagement  period  (August, 
1986;  Cornell,  1981;  Savage,  1986)  and  can  be  represented  by  a  high  order 
polynomial  (Cornell,  1981;  Boyd,  1987).  For  the  coordinate  system  as  shown 
Fig.  1,  the  following  gear  mesh  coupling  equations  are  established  by  equating 
forces  and  moments  (Choy,  1989).  For  the  kth  gear  of  the  system,  summing 
forces  in  the  X-direction  results  in: 
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Summing  forces  in  the  Y-direction  results  in: 
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Summing  moments  in  the  Z-direction  results  in: 
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where  n  is  the  number  of  gears  in  the  system. 

Using  the  modal  expansion  approach  (Choy,  1987,  1988a,  1989,  and  1990), 
the  motion  of  the  system  can  be  expressed  as: 
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where  m  is  the  number  of  modes  used  to  define  each  motion. 

Using  modal  expansion  and  the  orthogonality  conditions,  (Choy  1989,  1990) 
with  the  bearing  forces  due  to  the  base  motion  expressed  on  the  right  hand 
side  of  the  equations  the  modal  equations  of  motion  for  the  gear-rotor  system 
(Choy  1989)  can  be  written  ass 

X-Z  equation 
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which  can  also  be  expanded  into  modal  parameters  as  where 
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The  gear  mesh  forces  and  moments  can  also  be  expressed  in  the  modal  form,  for 
the  kth  gear  with  the  gear  location  at  the  1th  node,  ass 

m  n  - 
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where  k  is  the  gear-shaft-system  number,  j  is  the  mode  number,  and  1  is 
the  station  location  of  the  gear  mesh. 

A  set  of  modal  equations  of  motion  can  also  be  written  for  the  gearbox 
(Choy  1989,  1990)  as: 

N(M  +[Cbx]{M  *[^l  (M  *[Vr{FBx(t)}  -[$bxf{[Cxx][*](Al  +  [Cxy][*]W 

+  [Kxx][*](A}  +  [Kxy][$](B}}  *  0 

Similarly,  the  Y-equation  can  be  written  as 
E1]  {  Sb  }  +[Cby](M  +[Ab](Bb}  +[*byr(FBy(t)l  -[$byf{[Cyx][*](A}  +  [Cyy][$](6l  o 

+  [Kyx][*]{A)  +[Kyy][*HB}}  =  0 
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and  the  Z -equation  as 


W{flh}  +  [cbx]{^)  *  M{Db}  -  [*b,f  {*!«<*)  } 


The  procedure  involves  the  solution  of  the  coupled  modal  equations  of 
motion  between  the  (gear)  rotors  and  the  gearbox  structure.  The  coupling 
effects  of  the  gear  mesh  and  the  bearing  supports  are  also  expressed  in  modal 
coordinates  such  that  the  global  equations  are  solved  simultaneously  in  modal 
form.  A  set  of  initial  conditions  for  both  displacement  and  velocity  of  the 
global  system  are  calculated  from  the  steady  state  conditions  at  the  rotor¬ 
bearing  systems  and  zero  vibration  at  the  gearbox.  The  modal  accelerations  A, 
B,  A  ,  A.  B  ,  and  D  of  the  system  can  be  evaluated  (Eqs.  (11),  (13),  (15), 
and  (19)  to  (21)).  A  variable  time  stepping  integration  scheme  (the  Newmark- 
Beta  Method)  is  used  to  integrate  the  modal  acceleration  to  evaluate  the  modal 
velocities  and  displacements  at  the  next  time  step.  A  regular  time  interval 
of  200  points  per  shaft  revolution  is  used  in  this  study  except  for  a  refined 
region  of  smaller  steps  at  the  transition  from  single  to  multiple  tooth  con¬ 
tact.  The  modal  acceleration,  velocity,  and  displacement  calculated  from  the 
transient  integration  scheme  are  transformed  back  into  the  generalized  coor¬ 
dinates  (Bq.  (10)).  The  nonlinear  bearing  forces  and  gear  mesh  forces  can  be 
evaluated  from  the  velocity  and  displacement  differentials  between  the  rotors 
and  the  gearbox  structure. 


EXPERIMENTAL  STUDY 

The  gear  noise  rig  (Fig.  2)  was  used  to  measure  the  vibration,  dynamic 
load,  and  noise  of  a  geared  transmission.  The  rig  features  a  simple  gearbox 
(Fig.  3)  containing  a  pair  of  parallel  axis  gears  supported  by  rolling  element 
bearings.  A  150-kW  (200-hp)  variable-speed  electric  motor  powers  the  rig  at 
one  end,  and  an  eddy-current  dynamometer  applies  power-absorbing  torque  at  the 
other  end.  The  test  gear  parameters  are  shown  in  Table  I. 

Two  phases  of  experiments  were  performed  on  the  gearbox;  (1)  static  modal 
analysis  of  the  gearbox  and  (2)  dynamic  vibration  measurements  during  opera¬ 
tion.  Modal  parameters,  such  as  system  natural  frequenceis  and  their  corre¬ 
sponding  mode  shapes,  were  obtained  through  transfer  function  measurements 
using  a  two  channel  dynamic  signal  analyzer  and  modal  analysis  software.  For 
this  experiment,  116  nodal  points  (with  3°  of  freedom  each)  were  selected  on 
the  gearbox  housing.  Vibration  data  was  collected  from  accelerometers  placed 
on  the  gearbox  at  a  few  of  the  node  points  used  in  the  modal  survey.  Water¬ 
fall  plots  of  frequency  spectra  from  node  21  on  the  gearbox  top  and  from  node 
40  on  the  gearbox  end  are  presented  in  Figs.  4  and  5.  The  modal  frequencies 
(from  Figs.  6  and  7)  and  the  gear  mesh  frequency  are  shown  on  figures  4  and  5. 


DISCUSSION 

The  experimentally  obtained  modes,  shown  in  Figs.  6  and  7,  represent  the 
major  vibration  modes  of  the  gear  noise  rig  in  the  0  to  3  kHz  region. 

Although  these  modal  frequencies  are  only  a  small  part  of  the  total  modes  of 
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the  system,  they  represent  a  major  part  of  the  total  global  vibration  of  the 
system.  In  order  to  produce  a  compatible  analytical  simulation  of  the  test 
apparatus,  a  similar  set  of  modes  were  generated  using  a  finite  element  model 
of  the  gearbox  structure.  Out  of  a  total  of  25  modes  existing  in  the  analyti¬ 
cal  model  in  the  0  to  3  kHz  frequency  region,  the  eight  dominant  modes  were 
used  to  represent  the  gearbox  dynamic  characteristics.  These  analytically 
simulated  modes  are  shown  in  Figs.  8  and  9.  As  shown  in  Table  II,  the  natural 
frequencies  of  the  simulated  modes  are  within  5  percent  of  the  measured  modes. 
The  three-dimensional  analytical  mode  shapes  are  very  similar  to  the  experi¬ 
mental  modes  shapes  (Figs.  6  and  7).  The  correlation  in  the  results  between 
the  analytical  model  and  the  experimental  measurements  confirms  the  accuracy 
of  the  dynamic  representation  of  the  test  gearbox  using  only  a  limited  amount 
of  modes. 

In  the  experimental  vibration  spectra  of  Fig.  4  (from  the  gearbox  top), 
there  are  prominent  peaks  at  the  gear  mesh  frequency  on  the  traces  on  1500  rpm 
(700  Hz),  3750  rpm  (1750  Hz),  and  5500  rpm  (2567  Hz).  These  frequencies  are 
near  natural  frequencies  found  in  the  modal  survey  (658,  1762,  and  2536  Hz, 
respectively).  At  the  highest  speed  (5750  rpm),  the  higher  modes  (2536,  2722, 
and  2962  Hz)  are  excited  by  the  gear  mesh  vibration  and  its  sidebands. 

Similar  behavior  is  shown  in  the  vibration  spectra  from  the  side  of  the  gear¬ 
box  (Fig.  5).  In  both  of  these  figures,  the  amplitude  of  the  peaks  increases 
with  speed  as  the  mass  imbalance  force  increases  with  the  square  of  the  rota¬ 
tional  speed. 

The  vibration  response  of  the  gearbox  was  simulated  by  the  analytical 
method  presented  in  section  II.  The  gearbox  vibration  modes  were  calculated 
by  finite  element  methods.  The  gear-rotor  system  response  was  calculated  with 
a  discrete  rotor  model  using  the  matrix  transfer  technique.  The  gear-rotor 
system  was  coupled  to  the  gearbox  in  modal  coordinates  to  solve  for  the  tran¬ 
sient  vibration  response  of  the  global  system.  The  vibration  response  was 
then  transformed  into  the  frequency  domain  through  a  fast  Fourier  transforma¬ 
tion  routine. 

The  analytical  vibration  spectra  shown  in  Figs.  10  and  11  simulate  the 
experimental  results  of  Figs.  4  and  5.  Both  the  experimental  spectra  (Figs.  4 
and  5)  and  the  analytical  spectra  (Figs.  4  and  5)  show  excitation  at  the 
natural  frequencies  at  shaft  speeds  of  1500,  3500,  and  5500  rpm.  At  the 
higher  shaft  speed  (5500  rpm),  several  modes  are  excited  by  the  shaft  fre¬ 
quency  and  its  sidebands.  Small  differences  between  the  analytical  and 
experimental  spectra  may  be  due  to  (1)  small  errors  (<5  percent)  in  the  cal¬ 
culated  modes,  (2)  nonuniformaties  in  the  gearbox  not  present  in  the  modal, 

(3)  errors  in  measurement  of  rotational  speeds  and  vibration  amplitudes,  and 

(4)  limitations  in  modeling  the  bearings  as  a  radial  stiffness  element  only. 


CONCLUSIONS 

Analytical  and  experimental  studies  of  a  single  stage  gear  transmission 
system  were  performed.  Results  from  the  analytical  simulations  show  good 
correlation  with  experimentally  obtained  data  from  the  test  gearbox.  Some 
basic  conclusions  from  this  study  can  be  summarized  as  follows: 

1.  A  modal  synthesis  approach  can  be  used  to  simulate  the  dynamics  of  a 
gear  transmission  system. 
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2.  The  accurate  gearbox  model  can  be  developed  by  correlating  the  modal 
characteristics  from  experimental  study  with  those  from  analytical 
simulations . 

3.  The  proper  choice  of  modes  used  in  the  modal  synthesis  will  reduce 
the  number  of  modes  required  in  the  analysis,  without  sacrificing  accuracy  in 
the  numerical  solution. 
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TABLE  I. 


TEST  GEAR  PARAMETERS 


Gear  type  . 

Number  teeth  . 

Module,  mm  (diametrial  pitch  in."1)  .  .  . 

Face  width,  mm  (in.)  . 

Pressure  angle,  deg  . 

Theoretical  contact  ratio  . 

Driver  modification  amount,  mm  (in.)  .  . 
Driven  modification  amount,  mm  (in.)  .  . 

Driver  modification  start,  deg  . 

Driven  modification  start,  deg  . 

Tooth-root  radius,  mm  (in.)  . 

Gear  quality  . 

Nominal  (100  percent)  torque,  N-m(in.-lb) 


Standard  involute,  full-depth  tooth 

. 28 

.  3.174(8) 

.  6.35(0.25) 

. 20 

. 1.64 

.  0.023*0. 0009 ) 

.  0.025(0.0010) 

. 24 

. 24 

.  1.35(0.053) 

.  AGMA  class  13 

.  71.77(635.25) 


TABLE  II.  -  COMPARISON  OF  EXPERIMENTAL 
MEASURED  AND  ANALYTICAL  MODELED  NATURAL 


FREQUENCIES 


Dynamometer  - 


Figure  2. —  Picture  of  gear  noise  rig. 


Node  40- Y 
(on  back  side)  -x 
Node  21  -Z  \ 


Figure  3. —  Test  gear  box. 
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Modes  1  2  3  4  5  6  7  8 


Figure  4. —  Experimental  vibration  frequency  spectrum  at  node  21 
(note:  The  measured  modal  frequencies  are  shown  at  top  of 
figure) 


Modes  1  2  3  4  5  6  7  8 


Figure  5. —  Experimental  vibration  frequency  spectrum  at 
node  40  (note:  The  measured  modal  frequencies  are 
shown  at  top  of  figure). 
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Mode:  3  Mode:  4 

Frequency,  1709.77  Hz  Frequency,  1999.95  Hz 

Figure  6  —  Gear  box  experimental  mode  shapes. 
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Mode:  7  Mode:  8 

Frequency,  2722.16  Hz  Frequency,  2961  71  Hz 


Figure  7. —  Goar  box  oxporlmontal  modo  chapos 
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Mode:  1 

Frequency  -  658 


Mode:  2 

Frequency  -  1006 


Mode:  3 

Frequency  -  1762 


Mode  :  4 

Frequency  -  2051 


Figure  8  —  Gear  box  analytical  mode  shapes. 


Mode:  5 

Frequency  -  2336 


Frequency  -  2536 


y 


Mode:  7 

Frequency  *  2752 


Mode:  8 

Frequency  =  3012 


Figure  9. —  Gear  box  analytical  mode  shapes 
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V'b,all(  Amplitude,  rpm  Vibration  Amplitude,  rpm 


Mesh  frequency  -) 


Figure  10. —  Analytical  vibration  frequency  spectrum  at  node  21  (note:  The 
prediced  modal  frequencies  are  shown  at  top  of  figure). 


Modes  1 


Speed,  rpm 
5500  HLva-to-Aa.-' 


3  4  5  6  7  8 

i  i  mi 


4500 


3500 
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_ .  .  -A  A  A  ^ 


jIl 


2500 


1500 


. ...  )k. 

_ i 
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0  50  i  100 

Mesh  frequency-7 


150 


200 


250 


300 


350 


400 


Frequency.  kHz 

Figure  1 1 . —  Analytical  vibration  frequency  spectrum  at  node  40  (note:  The 
prediced  modal  frequencies  are  shown  at  top  of  figure). 
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